Atomistic insights into the binding of SARS-CoV-2 spike receptor binding domain with the human ACE2 receptor: The importance of residue 493.

SARS-CoV-2 is a coronavirus that has created a global pandemic. The virus contains a spike protein which has been shown to bind to the ACE2 receptor on the surface of human cells. Vaccines have been developed that recognize elements of the SARS-CoV-2 spike protein and they have been successful in preventing infection. Recently, the Omicron variant of the SARS-CoV-2 virus was reported and quickly became a variant of concern due to its transmissibility. This variant contained an unusually large number (32) of point mutations, of which 15 of those mutations are in the receptor binding domain of the spike protein. While several computational and experimental investigations comparing the binding of the Omicron and wild type RBD to the human ACE2 receptor have been conducted, many of these report contradictory findings. In order to assess the differential binding ability, we conducted 2 µs of classical molecular dynamics (cMD) simulation to estimate the binding affinities and behaviors. Based upon MM-GBSA binding affinity, per-residue energy decomposition analysis, center of mass distance measurements, ensemble clustering, pairwise residue decomposition and hydrogen bonding analysis, our results suggest that a single point mutation is responsible for the enhanced binding of the Omicron mutant relative to the WT. While the 15-point mutations in the receptor binding domain contribute positively and negatively to the affinity of the spike protein for the human ACE2 receptor, it is the point mutation Q493R that confers enhanced binding while the Q493K mutation results in similar binding. The MM-GBSA binding estimations over a 2 µs trajectory, suggest that the wild type binds to ACE2 with a value of -29.69 kcal/mol while the Q493K and Q493R Omicron mutants bind with energy values of -26.67 and -34.56 kcal/mol, respectively. These values are significantly different, given the error estimates associated with the MM-GBSA method. In general, while some mutations increase binding, more mutations diminish binding, leading to an overall similar picture of binding for Q493K and enhanced binding for Q493R.

Investigating Novel Thiazolyl-Indazole Derivatives as Scaffolds for SARS-Cov-2 MPro Inhibitors

COVID-19 is a global pandemic caused by infection with the SARS-CoV-2 virus. Remdesivir, a SARS-CoV-2 RNA polymerase inhibitor, is the only drug to have received widespread approval for treatment of COVID-19. The SARS-CoV-2 main protease enzyme (MPro), essential for viral replication and transcription, remains an active target in the search for new treatments. In this study, the ability of novel thiazolyl-indazole derivatives to inhibit MPro is evaluated. These compounds were synthesized via the heterocyclization of phenacyl bromide with (R)-carvone, (R)-pulegone and (R)-menthone thiosemicarbazones. The binding affinity and binding interactions of each compound were evaluated through Schrödinger Glide docking, AMBER molecular dynamics simulations, and MM-GBSA free energy estimation, and these results were compared with similar calculations of MPro binding various 5-mer substrates (VKLQA, VKLQS, VKLQG) and a previously identified MPro tight-binder X77. From these simulations, we can see that binding is driven by residue specific interactions such as π-stacking with His41, and S/π interactions with Met49 and Met165. The compounds were also experimentally evaluated in a MPro biochemical assay and the most potent compound containing a phenylthiazole moiety inhibited protease activity with an IC50 of 92.9 μM. This suggests that the phenylthiazole scaffold is a promising candidate for the development of future MPro inhibitors.